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ABSTRACT 

We study the formation and evolution of Brightest Cluster Galaxies starting from a 
z = 2 population of quiescent ellipticals and following them to z = 0. To this end, we 
use a suite of nine high-resolution dark matter-only simulations of galaxy clusters in 
a ACDM universe. We develop a scheme in which simulation particles are weighted 
to generate realistic and dynamically stable stellar density profiles at z = 2. Our ini- 
tial conditions assign a stellar mass to every identified dark halo as expected from 
abundance matching; assuming there exists a one-to-one relation between the visible 
properties of galaxies and their host haloes. We set the sizes of the luminous com- 
ponents according to the observed relations for z ~ 2 massive quiescent galaxies. We 
study the evolution of the mass-size relation, the fate of satellite galaxies and the 
mass aggregation of the cluster central. From 2 = 2, these galaxies grow on average 
in size by a factor 5 to 10 of and in mass by 2 to 3. The stellar mass growth rate 
of the simulated BCGs in our sample is of 1.9 in the range 0.3 < z < 1.0 consistent 
with observations, and of 1.5 in the range 0.0 < z < 0.3. Furthermore the satellite 
galaxies evolve to the present day mass-size relation by z = 0. Assuming passively 
evolving stellar populations, we present surface brightness profiles for our cluster cen- 
trals which resemble those observed for the cDs in similar mass clusters both at z = 
and at z = 1. This demonstrates that the ACDM cosmology does indeed predict mi- 
nor and major mergers to occur in galaxy clusters with the frequency and mass ratio 
distribution required to explain the observed growth in size of passive galaxies since 
z = 2. Our experiment shows that Brightest Cluster Galaxies can form through dissi- 
pationless mergers of quiescent massive z = 2 galaxies, without substantial additional 
star formation. 

Key words: galaxies: formation - galaxies: evolution - galaxies: clusters: general - 
galaxies: elliptical and lenticular, cD 



1 INTRODUCTION 

Brightest cluster galaxies (BCGs) form the massive end 
of the galaxy population. They are generally associated 
with old stellar populations, little star formation and large 
sizes jvon der Linden et al.ll2007l ; Ib ernar di 2009), except in 
strong cooling flows. In some clusters, BCGs are surrounded 
by a diffuse envelope of intracluster light. This additional 
light has been measu red in a number of nearby clusters 
(| Gonzalez et al.l 120051 ) as well as in stacks of BC Gs from 
the Sloan Digital Sky Survey (|Zibetti et al.l 120051 ) . There 
have been claims that the observed evo lution of BCGs 
(e.g. ICollins et all 120091 : IStott et all 120111 ) disagrees with 
the predictions of semi-analytic models of galaxy formation 



(IDe Lucia fc Blaizodl2007h . These studies selected BCGs in 
high redshift clusters and compared them to the central 
galaxies of present-day clusters of the same X-ray luminos- 
ity, observing little change in the sizes and stellar masses. 
However, recent results from Lidman et al. (2012) seem to 
indicate less tension between the models and observations. 
It is important to note that the current semi-analytic models 
of galaxy formation do not predict surface brightness profiles 
in a self-consistent manner. Thus, any direct comparison of 
sizes in the models and in real galaxies should be made with 
caution. 

Observations of z = 2 galaxies by iDaddi et all (2005) 
and iTruiillo et all l|2007t ) revealed the presence of a popula- 
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tion of massive quiescent galaxies with much smaller sizes 
than similar mass present-day elliptica ls. This result has 



been confirmed by several other groups (van Dokkum ct al 



2008; Ne wman et alj|2012f ). Minor mergers (jBezanson et al. 
2009 ; iNaab et al] 12009) have been proposed as the main 



mechanisms driving the recent size evolution of elliptical 
galaxies. Currently, it is unclear whether all su ch ga laxies 
will grow into present-day ellipticals. iBernardil (2009) pro- 
posed that some of these objects might be progenitors of 
tod ay's BCGs. 

IWhitd (|l97rj ) and lOstriker fc Hausmanl (|l977l 1 intro- 
duced galactic cannibalism as a possible mechanism to ex- 
plain the formation of cD galaxies: a central galaxy gradually 
swallows its companions as dynamical friction brings them 
to the cluster centre. These deposit their stellar material pri- 
marily on the outskirts of the larger galaxy, helping it grow 
in size and mass. This phenomenon seems to fit well within 
the ACDM cosmogony where structure grows hierarchically. 
For example, semi-analytic models of galaxy formation find 
that BCGs form in a two-phase process: an initial col- 
lapse with rapid cooling and star formation at high redshift 
is followed by later growth through mult iple dissipation- 
less mergers of pre-existing p rogenitors (|Khochfar fc Silkl 
120061 : |Pe Lucia fc Blaizoj 120071 ). A similar two-phase for- 
mation mechanism is also reported in hy drodynamical sim- 
ulations of massive elliptical formatio n llNaab et al] 120091: 
lOser et al.|[201ol : iFeldmann et alj|201ll ). Notably, lOser et all 
( 20121 ) show that compact massive ellipticals can grow onto 
the present-day mass-size relation through minor mergers. 
However, such simulations still produce galaxies that are too 
massive for the host haloes they inhabit so their quantitative 
app licability i s in so me doubt. 

iDubinskl (|l998l ) studied whether BCGs could form out 
of a population of spirals between redshift z — 3 and 
z = 0. His cosmologically motivated simulation showed 
that BCGs could potentially form through repeated dissi- 
pationless mergers of galaxies. This inspired later studies 
testing the collisionless merger hypothesis within cosmolog- 
ical simulations which provided new insights into the evo- 
lution of BCGs ijRuszkowski fc Sp ringe ]l2009l ) and into the 
origin of the intracluster light i Rudick et al.l 120061 ) . A per- 
sistent issue, however, is that the initial galaxies assumed 
appear imcons i stent with recent high redshift observations. 
iLaporte et ail (|2012l ) showed that re-scaling the luminous 
components of the galaxies to bring them into better agree- 
ment with these observations changed the size growth of 
BCGs from r oc M to r oc M 2 consistent with observational 
studies of massive quiescent galaxies at fixed number density 
ijvan Dokkum et al'llioiol : IPatel et al.ll2012l ) 

Here, we present a scheme which can weight particles 
in a high-resolution cosmological simulation of cluster for- 
mation to represent realistic stellar density profiles for the 
initial galaxies. Using this method we can assign to every 
halo of mass M at a certain redshift z%, a stellar mass m, 
according to an appropriate abundance matching relation , 

" Vale fc Ostrikeij | 2004| : 
Behroozi et all 12010 . 
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20131 ). This ensures a faithful represen- 
tation of the high redshift luminosity function at z = 2. 
We can then study the evolution of cluster galaxies from a 
population of galaxies consistent both with the luminosity 
function and with the mass-size relation observed at z — 2. 



We a i m to test the dissipationless merger hypoth esis (|Whitei 
1 19761 : lOstriker fc Hausmanlfl977l : |Pubinskilll99i ) . The ques- 
tions we want to address are the following: are BCGs special 
or simply the product of merging of normal cluster galaxies? 
Does star formation contribute substantially to the growth 
in mass of BCGs between 2 = 2 and z — 0? Do mergers 
explain why BCGs seem to lie off the present-day mass-size 
relation? Could BCGs have evolved from the observed high- 
redshift population of massive quiescent galaxies? 

Section 2 presents the simulations we use and our 
method of generating stellar density profiles. In section 3, 
we study the properties of the BCGs and compare them 
with observations. In section 4 we study the evolution of 
BCGs and how it compares to the population of cluster el- 
lipticals. We discuss the significance of our results in section 
5 and conclude in section 6. 



2 METHODS 
2.1 Simulations 

We use a set of nine dark- matter-only zoo m-in simulations o f 
galaxy clusters from the Phoenix project (|Gao et al] 120121 ) . 
These are named Ph-A to Ph-I. The haloes were initially 
select ed from the Millennium Simulation jSpringel et al] 
120051 ) and re-simulated with comoving softening length e = 
0.3/i -1 kpc and mass resolution m p ~ 4 — 10 x 10 6 fe~ 1 Mp ) . 
Details of the simulations are given in lGao et al] (2012). The 
subhaloes were ident ified with the structure finder Subfind 
jSpringel et al]|200ll ). We also compute the potential of ev- 
ery particle in each subhalo at redshift z — 2. For the rest 
of the discussion, our units are kpc and Mq for length and 
mass respectively. 



2.2 A Weighting scheme for cosmological dark 
matter simulations 



The weighting scheme pre sented here generalises that of 
iBullock fc Johnston! (|2005l ) for direct use in cosmological 
dark matter simulations of structure formation. Cold dark 
matter (just like stars in a galaxy) is collisionless and its 
distribution function (df) satistfies the collisionless Boltz- 
mann equation (CBE). Provided that the df of the dark 
matter in a halo is in a steady state, Jeans' theorem guar- 
antees the existence of a distribution function of the form 
/ = /(7i, /2,/3), where /i,2,3 are isolating integrals of the 
motion. In a triaxial potential a df of the form / — f(H), 
where H = \v 2 + $ is the Hamiltonian, can always be ger- 
enated. Generally, dark haloes are not in a steady state; 
their lives are continuously shaped by accretion events such 
mergers and infall. Howeve r, in their study of the e volution 
of the df of CDM haloes, iNataraian et aL I l|l997l) showed 
that between merger events haloes are in phases of "quasi- 
cquilibrium" within the virialised regions. Thus, provided 
we restrict ourselves to the virialised regions, Jeans' the- 
orem can be invoked to generate a stellar df of the form 
/. = /„(£). 

In order to generate a luminous stellar profile, we take 
each simulation particle of energy E to simultaneously rep- 
resent dark matter and stars in diferent amounts, through 
a weight function uj(E) = = ! * jyy 1 where N is 



N(E) 
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the differential energy distribution, g is the density of states 
and asterisks denote stellar quantities. 

We choose the s tellar distribution to be represented by 
the lHernauistl dl99Ch profile: 

(1) 



r (r + a) 

where a is the scale radius which is related to the half-mass 
radius through a = r e /(v2 + 1). The half-mass radius in 
projection is related to the half mass radius through R e = 
1.33r e 

We generate /*, using Eddington's formula: 

/.is) = -^-r r ^^.^^ + 



1 dv 
71 



(2) 

*=o 

where * = — $ + $ and £ = — £ + $ = W - u 2 /2 are the 
relative potential and total energies respectively. The poten- 
tial of an NFW profile is generally a reasonable description 
of the true potential and tends to zero in the limit r goes to 
infinity thus we set $o = 0. 

The density of states is given by: 

g(E) = (4tv) 2 f E r 2 ^j2{E-®(r))dr. (3) 
Jo 

Equations (2) and (3) are only valid for spherical sys- 
tems, however provided an educated guess, they can also be 
applied to triaxial ones. Indeed, authors in the past have 
already used the same formalism to study th e distribution 
function of dark haloes l|Nataraian et al]|l997l ) . To compute 

,2 

the term in Eddington's formula, we have to relate $ 
to v monotonically as <E>(r) is multivalued in a triaxial po- 
tential. For this, we simply approximate $(r) =< "l?(r) >. 
However, in the definition of the energy, we retain the actual 
value of the potential associated with every particle. This 
ensures the generation of a profile with equidensity surfaces 
which follow the potential. 

Figure 1 illustrates the kind of profile one can gener- 
ate within the z = 2 haloes of the Phoenix simulation. We 
have checked for the stability of our method by evolving a 
live dark halo for 150 Myr. We further looked for isolated 
galaxies in haloes which have not experienced any stellar 
mass growth between z = 2 and z = 0. We present such 
an example in Figure 2. This halo has grown in dark mat- 
ter mass by a factor of two but its stellar mass profile has 
remained unchanged over 10 Gyr. However, the half-light 
radius (defined as the radius of a sphere containing half the 
total stellar mass) has increased by 30 percent. Although 
this is not a substantial problem for studying the evolution 
of BCGs, where the change in radius is much larger and pri- 
marily caused by material deposited on the outside of the 
galaxy (|Laporte et alj|2012h it warns us that our experiment 
is close to the resolution limit for studying the size growth 
of galaxies which do not accrete much stellar mass. Consid- 
ering the large period of time over which this test has been 
carried out, we see that two-body relaxation processes are 
not the driving force behind the size growth we observe for 
the most massive galaxies. 

2.3 Initial Conditions 

For our initial conditions, we choose to represent a popula- 
tion of galaxies consistent with the observed stellar mass 
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Figure 1. Example pr 2 profile for a typical compact stellar pro- 
file generated using our weighting scheme. The dip is due to low 
particle numbers in the region below 1 kpc. Arrows mark the radii 
containing 50 and 95 percent of the stellar mass. 



function at z = 2, which we relate to ACDM haloes in 



the s imulations through recent AMR results (Most er et al 



Moster et al 



which are consistent with results from 
(2013) even at z = 2. We do not include the scatter in these 
relations, as the scatter for the mass-size relation is already 
large. We further require that the sizes of our model galax- 
ies are consistent with those observed for massive quiescent 
galaxies z — 2. For our pur poses, we choo s e the mass-size 



relation as parametrised bv I Williams et al and note 



Newman et al.l (|2012f) ( 



that it is very similar to that of [ 
the 3rd panel of their Figure 8). 

At z — 2, a few haloes are undergoing mergers, but be- 
cause these are a small fraction of the total we populate, 
we normally ignore them. However, if a particularly massive 
halo hosting a ~ 1O 11 M0 galaxy is undergoing a merger we 
generate its stellar distribution function at an earlier snap- 
shot. Depending on the state of the halo it is sometimes hard 
to generate stellar density profiles that perfectly match the 
original target, but as long as the galaxies are within the 
scatter of the observational relation we consider this good 
enough. For the rest of our discussion of the structural prop- 
erties of ellipticals in clusters we will focus on objects which 
have masses above 7 x 10 10 Mq, which still leaves us with a 
sample of 156 galaxies at z = 2. We follow, the subsequent 
evolution of galaxies by tracking the particle weights as the 
simulations evolve. 



3 STRUCTURAL PROPERTIES OF BCGS 

As a first test, it is interesting to check whether the total 
stellar mass of the final merger remnant agrees with expec- 
tations from the BCG luminosity function at z = 0. Fig- 
ure 3 shows the stellar-to-halo mass relation for the galax- 
ies at z = 2 and for cluster centrals and satellite galaxies 
at z = 0. As expected, many cluster galaxies have their 
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Figure 3. Stellar-to-halo mass relations at z=2 and z=0 (black 
and red lines respectively). The black stars represent the most 
massive progenitors of BCGs at z=2 and the squares represent 
all the haloes populated with stars in the initial conditions. The 
red triangles are the BCGs at z=0 and the red crosses represent 
cluster satellite galaxies. 



Figure 2. Top:pr 2 profile for a galaxy which has evolved from 
z = 2 (dashed line) to z = (solid line) without experiencing 
mergers or interactions with other haloes populated with stars. 
The host dark matter halo has grown in mass and half-mass radius 
by a factor of 2. This does not, however, translate in significant 
size growth of the stellar component. The arrows mark the radii 
enclosing 50 percent of the mass of the galaxies. Bottom: pr 2 
profiles for a BCG most massive progenitor (dashed black line), 
and final 2 = BCG (red line). We also show the distribution 
of the stars in the BCG most massive progenitor at z = (solid 
black line). The BCG grew by a factor of four in stellar mass and 
nine in size. The black and red arrows show the radii enclosing 50 
percent of the mass of the most massive BCG progenitor at z = 2 
and BCG at z = respectively. Clearly, the size growth observed 
for these galaxies is not due to two-body relaxation processes. 



dark matter haloes stripped, moving them horizontally in 
the m„ — Mhaio plane. Some of the cluster galaxies have a 
stellar mass deficit which can be explained by the lack of 
star formation in our experiment. These galaxies lie below 
the Moster et al. (2010) relation at z = 0. Some cluster 
satellite galaxies however grow to stellar masses in agree- 
ment with the relation at z — 0. Moreover, the final BCGs 
occupy a region that is in good agreement with expecta- 
tions of AMRs at z = 0. This success depends on the as- 
sumption we made about the progenitor galaxies at redshift 
z — 2 and on the hierarchical growth of the clusters accord- 
ing to the ACDM paradigm. At face value, this strongly 
supports the hypothesis that BCGs form from dissipation- 
less mergers of galaxies that were already in place at 2 = 2. 
This is in line with the galactic cannibalism picture origi- 



nally formulated by I White! (|l97rj| ) and lOstriker&T Hausmanl 
(1977): BCGs predominantly grow through the later assem- 
bly of already existing galaxies. This idea agree s with the 
more recent studiy of De Lucia fc Blaizoj (|2007l ) who find 
using their semi-analytic model that 80 percent of the stars 
ending up in BCGs are already formed at z = 3. With this in 
mind, we now ask whether any of the structural properties 
of our BCGs closely resemble those of known cluster central 
galaxies. 

3.1 Surface Brightness and Density Profiles 

Recent studies of ma s sive quiescent galaxies a t z — 2, 
(|Williams et al.l l20ld ; Ivan Dokkum et all l200ct ) adopt a 
Kroupa IMF and solar metallicity when interpreting their 
photometric data and find a stellar population age of 1 Gyr. 
Our experiment assumes that the stellar populations evolve 
passively (the collisionless merger hypothesis), so we deter- 
mine the ma ss-to-light r atios at later times using the values 
compu ted bv lMarastonl (|200 5,1. We note that lNewman et al.l 
il2012T) use Salpeter IMFs to derive their stellar masses, how- 
ever due to the large scatter in the mass-size relation this is 
a minor issue. In fact, the two sets of data agree at z = 2. 

In Figure 4, we present surface brightness profiles for 
BCGs at z = 0. This is done by taking 50 random pro- 
jections of the i ndividual galax i es. Su rface brightness pro- 
files taken from iGonzalez et al.l (|2005l) are overlaid in red. 
These authors measured the surface brightness profiles of 
24 nearby BCGs in the I-band. The clusters they used were 
in the redshift range 0.03 < z < 0.1. We find a reasonable 
match between the observed and simulated light distribu- 
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Figure 4. Surface brightness profiles for BCGs at z = derived 
assuming passive evolution between 2 = 2 and z = 0. Overplot- 
ted in red are surface brightness profiles for nearby galaxy clusters 
from Gonzalez (2005). The simulated surface brightness profiles 
have a sligthly steeper fall-off at large radii compared to the ob- 
served ones. 

tions both in shape and in normalisations. However, we still 
note that at large radii our simulated galaxies have system- 
atically slightly steeper fall-offs in their surface brightness 
profile. 

Some authors have claimed that BCGs are already well 
in place at high redshift z = 1 . and that they e volve little 
thereafter jCollins et alj|2009l ; Istott et alj|201ll ). We com- 
pare the surface brightness profiles o f some of our rela xed 
cluster BCGs to those of compiled bv lStott' et all (|201lh at 
z — 1.0 in Figure 5. These authors deduced a stellar popula- 
tion formation age at z — 3, similar to that is found for qui- 
escent z — 2 galaxies. The observations were performed in 
the HST F850LP band, but we compute our surface bright- 
ness profiles assuming a mass-to-light ratios in the SDSS z 
band which matches closely the F850LP band. We add the 
(1 + z) 4 surface dimming correction for direct comparison 
with their data. Given the uncertainties in the IMF and the 
different band zero point, the agreement between the two is 
still reasonable and encouraging. 

We also present the evolution in density profiles from 
z = 2 to z = in Figure 6, separating the components into 
in-situ and accreted components. As previous studies have 
shown the size evolution is predominantly drive n by adding 
stellar mass to the outskirts of the galaxies dNaab et al.l 
120091 ; lOser et alj|20l3 : lLaporte et aHl2012l ; iHilz et al.ll2012r i. 



4 EVOLUTION OF BCGS AND ELLIPTICALS 
IN CLUSTERS 

Recentlv. lLidman et alj (2012, hereafter L12), presented new 
results on the evolution of the stellar mass in BCGs using 




R [kpc] 

Figure 5. Simulated surface brightness profiles in the SDSS z- 
band for BCGs at z = 1 against projected radius. Overplotted 
in red are the observed surface brightness profiles in the F850LP 
HST band of for clusters of similar masses taken from the sample 
of Stott et al. (2011) 

samples of high, intermediate and low redshift clusters. We 
compare their data with ours in Figure 7 looking at cluster 
mass versus BCG stellar mass at three different redshifts. 
Our BCGs he within the scatter of their data, although gen- 
erally towards the massive end in stellar mass. However, we 
show this is also partly due to the way L12 derive their stel- 
lar masses. For some of our BCGs, we used de Vaucouleurs 
profile fits within an aperture of 30 kpc to derive stellar 
masses. For high Sersic indices (e.g. n — 8) up to half of the 
total stellar mass of the BCG can be missed. We indicate the 
shifts such systematics can cause for a few of our galaxies in 
Figure 7, showing that this might bring them in better agree- 
ment with the data. On the whole, L12 observe a stellar mass 
growth rate of 1.8 between 0.2 < z < 0.9 for their sample. 
This is consistent with that found for our simulated galax- 
ies: we observe a stellar mass-growth rate of 1.9, 1.5 and 2.6 
between 0.3 < 2 < 1.0, 0.0 < z < 0.3, and < z < 1 respec- 
tively. The stellar mass growt h since z < 1 report e d her e 
is consistent with results from |Pe Lucia fc Blaizotl l|2007t ); 
iTonini et all (|2012t ) 

We now turn to a comparison of the evolution in mass 
and size of the BCGs and that of the cluster population 
of ellipticals. We proceed by identifying all the z = 2 haloes 
which end up in the final clusters at z — (defined as the ma- 
jor FOF group) and track them over redshift using the trees 
constructed from the SUBFIND data at ~ 50 intermediate 
outputs using individual par ticle IDs to match subhaloes i n 
neighbouring outputs. As in Ruszkows ki fc Springell (2009), 
we compute three-dimensional spherical stellar mass profiles 
using logarithmic bins of widths A = 0.1 and measure their 
3D half-mass radii. This is shown in Figure 8. We visually 
inspected the individual density profiles of cluster galaxies 
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Figure 6. Density profiles of BCGs for the nine cluster simulations. We separate in-situ components at z = 2 (for the most massive 
progenitor) and z = (green and blue respectively) and accreted component at z = (magenta), the final BCG density profile is shown 
in black. Note that haloes B, H, I are unrelaxed at z=0 making their analysis cumbersome. 



checking their Hernquist profile fits given the parameters de- 
termined in the previous step, namely M and a. However, 
for BCGs, as their profile can sometimes be more extended 
(with Sersic profiles reaching n — 8 — 10) than a simple de 
Vaucouleurs law, these fits were often poor. Furthermore, 
we note that some clusters are out of equilibrium at various 
times which complicates the analysis of the central galaxy, 
this is the case for Ph-A-B-H-I at z = 1 and for the Ph-B-H-I 
clusters at z = 0. 

Comparing panels, we see that BCGs evolve more 
rapidly than other galaxies and that different BCGs evolve 
at different rates. While some of the BCGs experience rapid 
growth in stellar mass by z — 1 (e.g. Ph-A and Ph-H) , others 
see most evolution between z — 0.3 and z = (e.g. Ph-B, 
Ph-D Ph-G). There are also BCGs that show almost no sig- 



nificant evolution in stellar mass or size between z = 0.3 and 
z = (e.g. Ph-A, Ph-E). Furthermore, some BCGs still suf- 
fer major mergers between z — 0.3 and z — 0, this is the case 
for Ph-I, which doubles in st ellar mass. Such late phases of 
majo r merging are observed (|Liu et al.l 120091 ; iBrough et al.l 
l201ll) . We also see a growing mass gap between BCGs and 
other cluster ellipticals as we get closer to z — 0. 

Our result confirms that of iRuszkowski fc Spri ngcl 

(2009): BCGs tend to get off the mass-size relation of lo- 
cal ellipticals because of their much higher merger rate. 
Whether all BCGs should lie off the normal relation is 
harder to ascertain with a sample of only nine galaxy clus- 
ters. We find that at z = 0, our cluster-ellipticals and 
some our BCGs l ie on the mass-size relation derived by 
iHvde fc Bernardil (2009). For comparison we also show in 
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Figure 7. Cluster mass against BCG stellar mass at different 
redshifts. Top: High-z BCG sample (0.8 < z < 1.5) and the 
Phoenix BCGs at z = 1.0 (coloured stars). Middle: Intermediate- 
2 L12 BCGs (0.2 < z < 0.5) and Phoenix BCGs at Z = 0.3. 
Bottom: Low-z (z < 0.2) L12 BCGs and Phoenix BCGs. At z = 0, 
we show for some BCGs the stellar mass predicted by fitting a 
de Vaucouleurs profile within a 30 kpc aperture. For large indices 
(n=8 or n=10) this can lead to an underestimate of the true 
stellar mass of up to a factor of two. Models and observations are 
consistent at all z given the large scatter in the observations. 



Figure 7 the lShen et ail (|2003h relation which is shallower, 
a consequence of using Petrosian-based quantities which in- 
troduces a bias for objects of large Sersic index. 

It is interesting to also look at the amount of mass de- 
posited through mergers and disruption of satellites. We de- 
fine a merger event with the BCG when half of the stellar 
mass of a progenitor galaxy gets incorporated within the first 
subhalo of the FOF group (defined as the halo inhabited by 
the BCG). Anything below this threshold we called "diffuse" 
stellar accretion from stripped satellites. We list in Table 1 
the amount by which each type contributes to the mass ag- 
gregation for each BCGs. This table shows that for the vari- 
ous BCGs not all the mass accretion came from mergers but 
also some from stripped galaxies or galaxies in the process of 
merging. This amount of diffuse mass deposition, defined in 
Table 1 as the percentage of the final BCG mass diffusely ac- 



Run 


mo 


^acc 


^merger 


/diffuse 


M200 




1O 11 M 


1O 11 M 


10 11 M Q 




10 14 M Q 


Ph-A-2 


4.4 


8.6 


7.7 


0.07 


8.9 


Ph-B-2 


2.7 


7.4 


5.7 


0.16 


11.3 


Ph-C-2 


3.2 


7.5 


7.0 


0.05 


7.5 


Ph-D-2 


1.7 


3.8 


3.4 


0.07 


8.5 


Ph-E-2 


2.2 


9.0 


7.4 


0.14 


8.2 


Ph-F-2 


1.7 


5.6 


3.5 


0.20 


10.9 


Ph-G-2 


2.5 


5.8 


3.1 


0.32 


15.8 


Ph-H-2 


3.6 


3.5 


1.4 


0.30 


15.6 


Ph-I-2 


6.6 


16.4 


11.7 


0.20 


33.0 



Table 1. Table summarising initial mass of the most massive 
BCG progenitor, stellar mass accreted over the last 10 Gyr, stel- 
lar contributed by mergers, the fraction of the final stellar mass 
contributed by the disruption of satellites (not defined as merg- 
ers). Note that we define a merger to occur when more than 50 
percent of the stellar mass of a galaxy goes into the final BCG. 



Run 


1:1-2 


1:3-5 


1:10-100 


Ph-A-2 


2 


1 


4 


Ph-B-2 


1 


6 


1 


Ph-C-2 


1 


5 


5 


Ph-D-2 


1 


3 


5 


Ph-E-2 


3 





3 


Ph-F-2 


1 





2 


Ph-G-2 


1 





8 


Ph-H-2 





1 


2 


Ph-I-2 





6 


3 



Table 2. Table summarising the merger ratios for each BCG. 
Note that we define a merger when more than 50 percent of the 
stellar mass of a galaxy goes into the BCG. 



creted (i.e. fdi 



ffuse 



= (r 



~)/{mo + m a cc), where 



mo is the mass of the in-situ component of the most massive 
BCG progenitor), can vary between five percent (Ph-C) and 
thirty percent (Ph-G and Ph-H) of the total accreted stellar 
mass. Furthermore, these values of acc reted stellar mass ar e 
in agreement with the trends found in iMoster et aD (|2013h . 
We also present the number of mergers each BCGs experi- 
enced through their lifetime and list their ratios in Table 2. 
Although our definition of merger is somewhat arbitrary it 
clearly shows that BCGs go through a mixture and succes- 
sion of major (1:1 - 1:2) and minor merging events (1:3 - 
1:10 ) both of which are supported by observa tions of BCGs 
(e.g. lLiu et al.l2009l ; lE~dwards fc Pattonll2012T ). Furthermore, 
we see that the clusters for which we have identified fewer 
mergers have also a higher fraction of diffuse stellar accretion 
through stripped satellites: this is the case for clusters Ph- 
F, Ph-G and Ph-H. Although, there is a mixture of types of 
mergers, in the cases where minor merging has been predom- 
inantly acting (Ph-B, Ph-C, Ph-D, Ph-F, Ph-I), the numbers 
and ratios are in agreement wit h those found by the ad hoc 
simulations of iHilz et al] lj2013l ) to be needed to make the 
sizes of compact ellipticals grow to those of present day ones. 
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Figure 8. Mass-size evolution for cluster galaxies. Note that we use a 3D definition of half-mass radius for our simulated galaxies and we 
compare these to de-projected half-mass radii for observations. Black dots represent simulated ellipticals, stars depict the 2 = BCGs 
and their most massive progenitors in the higher redshift panels. Top left: z = 2 initial conditions overlaid with data from Newman et al. 
2012 (N12) in blue diamonds, van Dokkum et al. (2008) in red triangles and the relation from Williams et al. 2010. Top right: z = 1.0 
mass size relation with data from N12 in the redshift range 1.5 < z < 1.0. Bottom left: z = 0.3 mass-size relation with data from N12 in 
the redshift range 0.4 < z < 1.0. Bottom right: z = mass-size relation overlaid with the mass-size relations of Hyde & Bernardi (2009) 
and Shen et al. (2003) 



5 DISCUSSION 

We now turn to compare our work with prev i ous si mula- 
tion studies of BCGs. Recently, iMartizzi et al.l (2012) pre- 
sented hydrodynamical simulations of BCGs in which the 
central galaxy was brought to agree with the local AMRs 
through the inclusion of AGN feedback. However, the sur- 
face brightness profile of the final galaxy contains a stellar 
core of 10 kpc. Such large stellar cores are clearly absent 
in real BCGs, suggesting that the AGN feedback mech- 
anisms implemented in these simulations are too strong. 
While it seems that AGN feedback is the most promis- 



ing way of suppressi ng over-cooling and star formation at 
the centre of clusters JCrotpn et "aL 2006; Sii acki fc SpringeH 
120061 : IPuchwein et al.ll2010f ). its impact on the structure of 
the cluster central galaxy is probably negligible because too 
little gas is present in the star-dominated regions. Our work 
shows that BCGs can form through dissipationless mergers 
of galaxies between z — 2 and z — with properties closely 
resembling those seen in the local Universe (surface bright- 
ness profiles, stellar-to-dark matter ratios). 

Turning to the evolution of our by definition "quiescent" 
galaxies, we observe an evolutionary trend close to those 
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found in observations (e.g. iNewman et~a l. (2012)). However 
we caution that perhaps our experiment is not suited for 
such a direct comparison to observations. Indeed, the 2 = 2 
Universe is diverse, containing extended star-forming galax- 
ies as well as compa ct galaxies, both quies cent and star form- 
ing (e.g. Figure 2 in lSzomoru et all (|2012T l). In this paper we 
have assumed that all haloes at z = 2 contain compact qui- 
escent galaxies. On the other hand, we can affirm that the 
BCGs could well be descendants of quiescent massive ellipti- 
cals at z — 2. The best evidence for this in our experiment is 
that these galaxies started on the same mass-size relation as 
the bulk of galaxies. The high number of mergers they expe- 
rience turns them into present-day merger remnants which 
closely match the surface brightness profiles of known BCGs. 

One caveat of our experiment is that we omit the effect 
of the baryons on the total potential of galaxies at z = 2. 
We have considered the potential of the dark-matter only 
simulation to represent that of z — 2 galaxies. Including 
the self-gravity of the stars would in principle deepen the 
potential wells of the galaxies making tidal stripping of the 
inner-regions of haloes less efficient. However, the inclusion 
or omission of contraction of the dark matter haloes due to 
the presence of a stellar component in these simulations ap- 
pears to make relatively little difference to the final results 
on the evolution of the BCG (Ruszko wski fc Springel| [2009; 
lLaporte et afll2012h . Once more information is available on 
the internal structure of galaxies at z = 2, it will be neces- 
sary to consider more complex modelling including the stars 
explicitly. 



6 CONCLUSIONS 

We have studied the late formation and assembly of BCGs 
in the ACDM cosmology starting from a z — 2 a popula- 
tion of "quiescent" galaxies resembling those observed. We 
did so by representing the stellar components of galaxies in- 
side their dark matter haloes through a distribution function 
weighting scheme. In contrast with previous studies, our ex- 
periment simultaneously represents the luminosity function 
at z = 2 as derived from abundance matching relations and 
the observed mass-size relation at z — 2. Previous studies N- 
body studies of this have asssumed substantially too much 
stellar mass for a given h alo mass (|Ruszkowski fc Springell 
120091 ; iRudick et al.ll201lh . thus adopting initial conditions 
inconsistent with modern observations and a ACDM cos- 
mology. 

Under the assumptions of our experiment, we predict 
present-day BCGs with stellar masses in good agreement 
with those inferred from the z — stellar-to-halo mass 
relation. Moreover, the surface brightness profiles of our 
simulated BCGs match local observations of known BCGs 
both in shape and in normalisation. This suggests that most 
BCGs have evolved passively from z = 2 to z = 0, forming 
out of pre-existing a population of compact galaxies with 
very little star formation after z = 2. 

The large masses and sizes of cD galaxies reflect the 
hierarchical growth of structure in ACDM and their special 
location at the centres of galaxy clusters. Our experiment 
also reproduces su rface brightness p rofiles of some of the 
five z ~ 1 BCGs in lStott et al.l (|20 1 if ) suggesting that these 
do not conflict with ACDM expectations. In fact, we show 



that our results agree well with the recent study of Lidman 
et al. (2012). Our simulated galaxies suggest that estimates 
based on fluxes within 30 kpc apertures may in some cases 
substantially underestimate the stellar masses of BCGs. 

The Phoenix project offers an important window to 
study various aspects of cluster dynamics and the evolu- 
tion of galaxies. Our presented scheme could be used in the 
future with the aim to address questions about colour and 
metallicity structure in BCGs using newly developed semi- 
analytic models (Yates et al. in prep). 
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